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ABSTRACT 



The high-frequency quasi-periodic oscillations (HP QPOs) that appear in the X-ray fluxes of low-mass X-ray binaries remain an 
unexplained phenomenon. Among other ideas, it has been suggested that a non-linear resonance between two oscillation modes in 
an accretion disc orbiting either a black hole or a neutron star plays a role in exciting the observed modulation. Several possible 
resonances have been discussed. A particular model assumes resonances in which the disc-oscillation modes have the eigenfrequen- 
cies equal to the radial and vertical epicyclic frequencies of geodesic orbital motion. This model has been discussed for black hole 
microquasar sources as well as for a group of neutron star sources. Assuming several neutron (strange) star equations of state and 
Hartle-Thorne geometry of rotating stars, we briefly compare the frequencies expected from the model to those observed. Our com- 
parison implies that the inferred neutron star radius Rf^g is larger than the related radius of the marginally stable circular orbit r,^s for 
nuclear matter equations of state and spin frequencies up to 800Hz. For the same range of spin and a strange star (MIT) equation of 
state, the inferrred radius is Rf^s ~ ''ms- The "Paczyhski modulation" mechanism considered within the model requires that Rf^s < ''ms- 
However, we find this condition to be fulfilled only for the strange matter equation of state, masses below IMq, and spin frequencies 
above 800Hz. This result most likely falsifies the postulation of the neutron star 3:2 resonant eigenfrequencies being equal to the 
frequencies of geodesic radial and vertical epicyclic modes. We suggest that the 3:2 epicyclic modes could stay among the possible 
choices only if a fairly non-geodesic accretion flow is assumed, or if a different modulation mechanism operates. 
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1. Introduction: HF QPOs and desire for 
strong-gravity 

Galactic low mass X-ray binaries (LMXBs) display quasi- 
periodic oscillations (QPOs) in their observed X-ray fluxes (i.e., 
peaks in the X-ray power density spectra). Characteristic fre- 
quencies of these QPOs range from ~10"^Hz to -lO^'Hz. Of 
particular interest are the so-called high-frequency (HF) QPOs 
with frequencies typically in the range 50- 1300 Hz, which is 
roughly of the same order as the range of frequencies character- 
istic of orbital motion close to a low mass compact object. We 
briefly recall that there is a crucial difference between HF QPOs 
observed in black hole (BH) and neutron star (NS) systems. In 
BH systems, the HF QPO peaks are commonly detected at con- 
stant (or nearly constant) frequencies that are characteristic of a 
given source. When two or more QPO frequencies are detected, 
they usually come in small-number ratios, typically in a 3 : 2 ra- 
tio (Abramowicz & Kluzniak, 2001; Kluzniak & Abramowicz, 
2001; McClintock & Remiflard, 2004; Torok et al., 2005). For 
NS sources, on the other hand, HF (or kHz) QPOs often appear 
as twin QPOs. These features, on which we focus here, consist of 
two simultaneously observed peaks with distinct actual frequen- 
cies that substantially change over time. The two peaks forming 
twin QPO are then referred to as the lower and upper QPO in 
agreement with the inequality in their frequencies. 

The amplitudes of twin QPOs in NS sources are typically 
much stronger and their coherence times much higher than those 
in BH sources (e.g. McClintock & Remillard, 2004; Barret et al., 
2005a,b; Barret et al., 2006; Mendez, 2006). It is however in- 
teresting that most of the twin QPOs with high statistical sig- 



nificance have been detected at lower QPO frequencies 600- 
700Hz and upper QPO frequencies 900- 1200Hz. Because of 
this the twin QPO frequency ratio clusters mostly around x; 3:2 
value posing thus some analogy to BH case (see Abramowicz et 
al., 2003a; Belloni et al., 2007; Torok et al., 2008a; Torok et al., 
2008b, c; Boutelier et al., 2009, for details and a related discus- 
sion). In several NS sources, the difference in the amplitudes 
of the two peaks changes sign as their frequency ratio passes 
through the (same) 3:2 value (Torok, 2009). A detailed review 
on the other similarities and differences in the HF QPOs features 
can be found in van der Klis (2006). 



1.1. HF QPO interpretation 

There is strong evidence supporting the origin of the twin QPOs 
inside 100 gravitational radii, rg = GMc"^, around the accreting 
compact objects (e.g., van der Klis, 2006). At present, there is 
no commonly accepted QPO theory. It is even unclear whether 
this theory could involve the same phenomena for both BH and 
NS sources. Several models have been proposed to explain the 
HF QPOs, most of which involve orbital motion in the inner re- 
gions of an accretion disc. When describing the orbital motion, 
the Newtonian approach necessarily fails close to the compact 
object. Two of the most striking differences arise from the rel- 
evant general relativistic description: Einstein's strong gravity 
cancels the equality between the Keplerian and epicyclic fre- 
quencies, and (due to the existence of the marginally stable cir- 
cular orbit Tnjs) it applies a limit to the maximal allowed orbital 
frequency. Several effects such as the relativistic precessions of 
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orbits then pop up in the inner accretion region. Finding a proper 
QPO model may thus help us to test the strong field regime pre- 
dictions of general relativity and, in the case of NS sources, also 
the models of highly dense matter (see van der Klis, 2006; Lamb 
& Boutloukos, 2007, for a review). 



1.2. Non-linear resonances between "geodesic and 
non-geodesic" disc-osciiiations 

Numerous explanations of the observed lower and upper HF 
modulation of the X-ray flux have been proposed while hypo- 
thetical resonances between the two QPO oscillatory modes are 
often assumed. Specific ideas considering non-linear resonances 
between disc-oscillation modes have been introduced and ex- 
tensively investigated by Abramowicz, Kluzniak and collabora- 
tors (Kluzniak & Abramowicz, 2001; Abramowicz & Kluzniak, 
2001; Abramowicz et al., 2003a,b; Rebusco, 2004; Torok et al., 
2005; Horak, 2008; Stuchlik et al., 2008; Horak et al., 2009, and 
others; see also Aliev & Galtsov, 1981 and Aliev, 2007). These 
ideas have been widely discussed and adapted into numerous in- 
dividual disc-oscillation models. 

The subject of disc oscillations and their propagation has 
been extensively studied analytically for thin disc (i.e., nearly 
geodesic, radiatively efficient) configurations (Okazaki et al., 
1987; Kato et al., 1998; Wagoner, 1999; Wagoner et al., 2001; 
Silbergleit et al., 2001; Ortega-Rodriguez et al., 2002; Wagoner, 
2008). The derived results have been compared to those for 
"thick" (radiatively ineflicient, slim-disc or toroidal) config- 
urations foe which both analytical (Blaes, 1985; Sramkova, 
2005; Abramowicz et al., 2006; Blaes et al., 2006, 2007; 
Straub & Sramkova, 2009) and numerical (Rezzolla et al., 
2003a,b; Rezzolla, 2004; Montero et al., 2004; Zanotti et al., 
2005; Sramkova et al., 2007) studies have been performed. 
Several consequences of disc-oscillation QPO models have been 
sketched, some having direct relevance to non-linear resonance 
hypotheses. In particular, it has been found that, due to pressure 
effects, the values of the frequencies at radii fixed by a certain 
frequency ratio condition can differ between the geodesic and 
fairly non-geodesic flow of factors such as 15% (Blaes et al., 
2007). 

1.3. Aims and scope of this paper 

Gondek-Rosiriska & Kluzniak (2002) suggested that the reso- 
nance theory of kHz QPOs can help us to discrimine between 
quark (strange matter) stars and neutron stars. In the spirit of this 
suggestion, we examine a particular, often quoted "3:2 epicyclic 
resonance model" (or rather a class of these models). The paper 
is arranged as follows. 

In Sect. 2, we briefly highlight some important aspects of 
non-linear resonance models specific to neutron stars and a 3:2 
epicyclic resonance model. In Sect. 3, we compare the model 
to the HF QPO observations of a group of NS sources display- 
ing the 3:2 ratio. The restrictions to the mass and radius implied 
by the equations of state for non-rotating NS are included. In 
Sect. 4, we explore the corrections required for NS rotation and 
again consider the equations of state. In Sect. 5, we assign some 
consequences and discuss possible falsification of the model, 
whereas the nearly geodesic and fairly non-geodesic cases are 
considered separately. 

Throughout the paper, we use the standard notation where 
Vi, Vu represent the observed lower and upper QPO frequencies, 
while Vjf , Vr, vg represent the Keplerian, radial epicyclic and ver- 



tical epicychc frequencies for the considered spacetime and its 
parameters. 



2. Resonances in discs around neutron stars 

Miscellaneous variations in the non-linear, disc-oscillation res- 
onances have been discussed in the past (see, e.g., Abramowicz 
& Kluzniak, 2001, 2004; Torok, 2005a; Horak & Karas, 2006). 
While the basic approaches have been common to both black- 
hole and neutron-star models, several differences between the 
two classes of sources have been considered. In particular, it has 
been suggested that, in a turbulent NS accretion flow, the res- 
onant eigenfrequencies are not fixed (e.g., when oscillations of 
a tori changing its position are assumed; Zanotti et al., 2003; 
Rubio-Herrera & Lee, 2005; Abramowicz et al., 2006; Torok 
et al., 2007; Kluzniak et al., 2007), or that the resonant cor- 
rections to eigenfrequencies reach high values (Abramowicz et 
al., 2003b, 2005a,b). Both possibilities are taken into account in 
Sect. 3. 



2.1. Modulation 

One more important difference between the two aforementioned 
source types concerns the QPO modulation mechanism (Bursa et 
al., 2004; Horak, 2005a; Abramowicz et al., 2007; Bursa, 2008). 
In the black hole case, the weak modulation is assumed to be 
primarily connected to radiation of the oscillating disc and the 
related relativistic lensing, light-bending, and Doppler effects. 
In the neutron star case, the expected modulation is connected 
to the flux emitted from a hot spot on the NS surface causing a 
strong QPO amplitude. 

We briefly describe the "Paczyriski-modulation" mechanism 
(Paczyhski, 1987), which was investigated by Horak (2005a) 
and Abramowicz et al. (2007). The schematic Fig. 1 displays 
the considered situation. The expected mass flow is described 
by the Bernoulli equation, while surfaces of constant enthalpy, 
pressure, and density coincide with surfaces of constant effec- 
tive potential 1A{r,z) - constant (Abramowicz, 1971). The disc 
equilibrium can exist if the disc surface corresponds to one of 
the equipotentials inside the so-called Roche lobe (region indi- 
cated by the yellow colour). No equilibrium is possible in the 
region of r < ^n. For a given accretion rate the dynamical mass 
loss occures when the fluid distribution overflows the surface of 
the disc for IAq — llir-i^). When the accretion disc oscillates, 
it slightly changes its position with respect to the equipotencial 
surfaces. At a particular location corresponding to the crossing 
of the equipotentials, the so-called cusp, even a small displace- 
ment of the disc causes a large change in the accretion rate. The 
change in the accretion rate is then nearly instantly reflected by 
the hot-spot temperature leading to an enhanced X-ray emission 
(Paczyhski, 1987; Horak, 2005a; Abramowicz et al., 2007). 

The existence of the surface I/q above the neutron star is 
crucial to the model. Therefore, as a necessary condition for its 
applicability, it is required that 



RNs/rms < 1, 



(1) 



where /?ns denotes the neutron (strange) star radius ("accretion 
gap paradigm", Kluzniak & Wagoner, 1985; Kluzniak et al., 
1990). We note that this is a necessary but insufficient condition, 
since the inner radius ^n is located between the marginally stable 
and marginally bound circular orbit (Kozlowski et al., 1978). 
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Fig. 1. Mass-flow leaving the disc and crossing the relativistic accretion 
gap (after Abramowicz et al, 2007). Top: Keplerian angular momentum 
versus the angular momentum in the flow. Bottom: The equipotential 
surfaces and the distribution of fluid in a meridional cross-section of 
the disc-configuration. The yellow area denotes the fluid in the disc, 
while the orange area corresponds to the overflow modulated by the 
oscillations. Enhanced luminosity arises as the flow enters the boundary 
layer (light-blue colour). 

2.2. Epicyclic resonance 

A particular example of the non-linear resonance between disc- 
oscillation modes is represented by the concept of the "3:2 
epicyclic internal resonance". This hypothesis is widely dis- 
cussed (e.g., Abramowicz et al., 2002; Kluzniak & Abramowicz, 
2002, 2005; Horak, 2004, 2005b; van der Klis, 2005; Torok & 
Stuchlili, 2005; Vio et al., 2006; Rebusco, 2008; Reynolds & 
Miller, 2009, among the other references in this paper). It as- 
sumes that the resonant modes have eigenfrequencies equal to 
radial and vertical epicyclic frequency of geodesic orbital mo- 
tion given by 



Vr('"3:2), V^, = Ve{r3;2), 



(2) 



associated with the orbital radius r^a, where vg/vi- = 3/2. 
We emphasize that models consider oscillations of fluid con- 
figurations rather than test particle motion (see, e.g., Kluzniak, 
2008, for some details and related references). In the following 
sections we consider the 3:2 epicyclic model and whether the 
Paczyiiski modulation mechanism may be at work. 

3. NS mass and radius implied by the 3:2 epicyclic 
resonant model 

In resonance models of BH QPOs, the observed constant fre- 
quencies are expected to coincide with the resonant eigenfre- 
quencies. Assuming a particular resonance, one may then relate 



the black hole spin or mass to the observed frequencies. This 
procedure was followed by Abramowicz & Kluzniak (2001) and 
later by Torok et al. (2005) and Torok (2005b) for various reso- 
nances and sets of sources. In principle, similar calculations can 
also be made for resonance models of NS QPOs. For neutron 
stars, the observed frequencies, however, change over time and, 
moreover, monotonic positive frequency correlations are similar, 
but specific to the individual sources. Within the framework of 
the resonance models we can consider two distinct simplifica- 
tions to the observed frequency correlations when inferring the 
neutron star mass: 

a) The observed frequencies are roughly equal to the resonant 
eigenfrequencies and the observed frequency correlation fol- 
lows from the changes in eigenfrequencies 



Vl = v!,'(r3:2 + Ar), v, = v^(r3:2 + Ar), 
implying for the 3:2 epicyclic model that 

Vl = Vi-(r3:2 -H Ar), V„ = Vfl(r3:2 + Ar). 



(3) 



(4) 



b) The eigenfrequencies are constant and the observed correla- 
tion is caused by the resonant corrections 



implying for the 3:2 epicyclic model that 

Vl = Vi-(r3:2) H- Av^, V„ = Vgiria) + Av„. 



(5) 



(6) 



We note that for a) the resonance plays rather a secondary role in 
producing QPOs, while for b) it represents their generic mecha- 
nism. 

3.1. Mass 

In this section, we neglect the effects of neutron star spin and as- 
sume the 3:2 epicyclic resonance model in the Schwarzschild 
spacetime.' Introducing a relativistic factor T = c^ I(2jiGM), 
Eq. (4) reads 



v^. = x'^''^ T, x = r/M, 



ve 



Vt = Vr S Vj. -\/l , 

X 



implying that 



-(^r- 



(7) 



(8) 



It has been previously discussed in terms of a correlation be- 
tween the QPO frequency (v^ or Vy) and frequency difference 
Av - v,j - Vl that the correlation given by Eq. (8) clearly dis- 
agrees with the observations of NS sources (e.g., Belloni et al., 
2005). The 3:2 epicyclic resonance model that is fully based on 
Eq. (3) is therefore excluded. Hence, in the following we focus 
on the option represented by Eq. (5). 

The relation of Eq. (5) to the observation of several NS 
sources was considered by Abramowicz et al. (2005a,b). They 
assumed that the corrections Av in Eq. (6) vanish when the ob- 
served frequency ratio Vu/vl reaches the 3/2 value. They sug- 
gested that the resonant eigenfrequencies [v^, vj'] in a group of 



' We consider here this standard spacetime description for non- 
rotating neutron stars, although some alternatives have been discussed 
in a similar context (see Kotrlova et al., 2008; Stuchlik & Kotrlova, 
2009). 
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twelve NS sources are roughly equal to [600Hz, 900Hz]. For the 
3:2 epicyclic model we then find that 



2.5 



.,3:2 



= 600Hz, 



,,3:2 



= 900Hz, 



(9) 



which, from terms given in equation (7), implies that the relevant 
mass must be around M - IMq (as first noticed by Bursa 2004 
unpublished). 

3.2. Radius 

Modelling of NS equations of state (EoS) have been exten- 
sively developed by numerous published methods and codes 
(see, Lattimer & Prakash, 2001 and Lattimer & Prakash, 2007 
for a review). Here we calculate NS radii following the approach 
of Hartle (1967), Hartle & Thorne (1968), Chandrasekhai- & 
Miller (1974), and Miller (1977). In Fig. 2, we plot the mass- 
radius relations for several EoS. 

Skyrme represents nine different EoS (namely SkT5, SkO', 
SkO, SLy4, Gs, SkI2, SkI5, SGI, and SV) given by the dif- 
ferent parameterizations of the effective Skyrme potentials (see 
Rikovska Stone et al., 2003, and references therein). DBHF rep- 
resents four different parameterizations, chosen to describe mat- 
ter in the framework of Dirac-Brueckner-Hartree-Fock theory. 
In particular, we choose the parameterizations labeled HA, HB, 
LA, and MA in Kotulic Bunta & Gmuca (2003) used by Urbanec 
et al. (2010) to describe the properties of static neutron stars. 
The EoS labeled APR has often been used. We chose the model 
labeled A\% + 6v + UIX* in the original paper (Akmal et al., 
1998). Remaining pure neutron-star equations of state are FPS 
(Pandharipande & Ravenhall, 1989) and BBB2 (Baldo et al., 
1997). The model labeled GLENDNH3 also includes hyperons 
(Glendenning, 1985). 

The MIT model represents strange stars calculated using the 
so-called MIT bag model (Chodos et al., 1974), where we used 
the standard values B = 10''*g.cm"'' for the bag constant and 
ffc =0.15 for strong interaction coupling constant. 

From Fig. 2 we can see that for M ~ IMq the NS radii are 
in all cases above r^s- Thus, assuming the Schwarzschild metric 
the condition presented in Eq. (1) is not fulfilled for the X-ray 
modulation given by the 3:2 epicyclic model. 

4. Effects related to NS spin 

We have restricted attention to the implications of Eq. (9) for 
non-rotating NS. The spin of the astrophysical compact objects 
and related oblateness however introduce some modifications of 
the Schwarzschild spacetime geometry. Without the inclusion of 
magnetic field effects, it has been found that the rotating space- 
times induced by most of the up-to-date neutron star equations 
of state (EoS) are well approximated with the solution of Hartle 
& Thorne (1968) (see Berti et al., 2005, for details). We use this 
solution (in next HT) to discuss the spin corrections to the above 
results. 

The HT solution reflects three parameters, the neutron star 
mass M, angular momentum J, and quadrupole moment Q. We 
note that the Kerr geometry represents the "limit" to the HT ge- 
ometry for q s QM/J^ -^ 1 up to the second order in /. The for- 
mulae for Keplerian and epicyclic frequencies in the HT space- 
time were derived by Abramowicz et al. (2003c). We applied 
these formulae to solve Eq. (9). Figure 3 displays the resulting 
surface colour-scaled in terms of M/Mq, j = cJ/GM^, and q. 
We can see that for low values of q and any j the implied M in- 
creases with increasing j, while exactly the opposite dependence 
M{j) occurs for high values of q and j > 0.2. 



1.5 



^ 
^ 



0.5 



1 1 1 1 


1 


DBHF ^^^^^^■ 

" APR ^^^^^^^ 

BBB2 XV/V^ 

GLEN - \Y \\M; 

MIT ^^/tTUV 


L 


1 1 1 1 


1 



10 
^Ns [km] 



12 



14 



Fig. 2. Mass-radius relations for several EoS assuming a non-rotating 
star The shadow area indicates the region with NS radii higher than 
the radius of the marginally stable circular orbit (no accretion gap). The 
mass M = I ± O.IMq is denoted by the dashed and dotted horizontal 
lines. 



4. 1 . EoS and radii 

For a given EoS, the parameter q decreases with increasing 
MjM^sj.- In more detail, it is usually q ~ \Q for IMq, while 
q € (1.5, 3) for the maximal allowed mass (e.g. Torok et al., 
2010, Fig. 3 in their paper). Since the non-rotating mass inferred 
from the model is about IMq, one can expect that realistic NS 
configurations will be related to M-j solutions associated with 
high q. These are denoted in Fig. 3 by colours of the yellow-red 
spectrum. 

We checked this expectation using the same set of EoS as in 
Sect. 3.2. We calculated the configurations for each EoS cover- 
ing the range of the central density p^ implying that M e (x; 
0.5Mo, M,nax) and the spin frequency Vs e (0, v,nax), using 
thousand bins in each of both independent quantities. The mass 
Mn,ax is the maximal mass allowed for a given EoS and spin fre- 
quency Vs. The frequency v^ax is the maximal frequency of a 
given neutron star and is equal to the Keplerian frequency at the 
surface of the neutron star at the equator, corresponding to the 
so-called mass-shedding limit. In this way, we obtained a group 
of 15 X 1000^ a; 10^ configurations. From these we retained only 
those fulfilling the condition in Eq. (9) for the epicyclic frequen- 
cies (v = V (M, j, q), Abramowicz et al., 2003c) extended to 



2/3v, 



,3:2 






(580Hz, 680Hz). 



(10) 



We note that this range of considered eigenfrequencies is 
based roughly on the range of the observed 3:2 frequencies 
(Abramowicz et al., 2005a,b). The combinations of mass and 
angular momentum selected in this way are displayed in Fig. 4a. 
Inspecting the figure, one can see that the mass decreases with 
increasing j above j ~ 0.3. 

Figure 4b indicates the ratio Rus/fms for the selected config- 
urations (shown in Fig. 4a). The modulation-condition presented 
in Eq. (1) appears only to be fulfilled for MIT-EoS and high spins 
above 800Hz. 

5. Discussion and conclusions 

The neutron star masses inferred for the 3:2 epicyclic resonance 
model by the considered EoS (Fig. 4a) are very low compared to 
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a) 



Fig. 3. Solution of the 3:2 frequency equation (9) projected onto M - j 
plane and colour-scaled in terms of q. 



the "canonical" value of IAMq. For the non-rotating case, the 
implied NS configurations are in addition insufliciently compact 
to fulfill the modulation condition in Eq. (1). We find that this 
condition is satisfied only for high spin values, above 800Hz, 
and strange matter EoS (MIT) (see the shaded region in Fig. 4b). 

Searching through the region, we find the highest mass satis- 
fying Eq. (1) to be M = 0.97 Mq. The related NS spin is 960Hz. 
This mass and spin correspond to v^'^ - 580Hz. For higher fre- 
quencies v^'--, the required mass is even lower. For v^'-^ - 630Hz, 
it is M = O.85M0, whereas the related NS spin is 900Hz. 

For compact objects in the NS kHz QPO sources, there are 
at present no clear QPO independent mass estimates. In contrast, 
there is convincing evidence of the spin of several sources from 
the X-ray burst measurements (see, e.g, Strohmayer & Bildsten, 
2006). In the group of sources discussed by Abramowicz et al. 
(2005a,b) considered in this paper, there are several that have 
spins in the range ~250-650Hz. The NS parameters implied 
by the 3:2 epicyclic model therefore include not only very low 
masses, but also spins excluded by the QPO independent meth- 
ods. 

The results obtained that falsify the epicyclic hypothesis are 
doubtless as far as: 

i) The Paczynski modulation mechanism is involved, implying 
that the inequality R^s < '"ms is vaUd. 

ii) The eigenfrequencies considered within the model are equal 
to (nearly) geodesic frequencies. 

The amplitudes of NS twin peak QPOs are often far higher 
than the BH amplitudes. The lensing effects are then insuffi- 
cient for the observed modulation (see e.g.. Bursa et al., 2004; 
Schnittman & Rezzolla, 2006). Nevertheless, there may be 
mechanisms other than Paczynski modulation (that do not nec- 
cessarily require the /?ns < /"ms condition). One example^ may 
be an equilibrium torus without a cusp that oscillates with a high 
oscillation amplitude. In this case, accretion onto the neutron star 
may yet occur due to the overflow of the critical equipotential. 
The largest configurations fullfilling the 3:2 frequency condition 
for a non-rotating NS have the radius 1 .8 x6M = 10. 8M, which 
equals the radius of the 3:2 resonant orbit for j - 0. Since all 



^ The authors thank L. Rezzolla for suggesting this possibility during 
a discussion at the Relativistic Whirlwind conference in Trieste (June 
2010). 
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Fig. 4. a) NS configurations fulfilling the 3:2 frequency condition (10). 
b) Related relationships between NS spin and radii evaluated in terms 
of ISCO radii r„,s(M, j, q). Only subset of configurations obeying the 
condition (10) are depicted in the figures for clearness. Lines tending to 
appear on both figures correspond to configurations with same central 
parameters and different rotational spin. 



the considered EoS infer a NS radius below the resonant radius, 
the model should, in principle, work for this hypothetical mech- 
anism. In addition, interference between a terminating disc and 
spinning NS surface near the 3:2 resonant orbit could represent 
a powerful excitation mechanism (see also Lamb & Coleman, 
2003). Nevertheless, some difficulties are apparent. In particu- 
lar, the (unexplored) mechanism should be adjusted to the flow 
coming from the binary companion. The twin QPO peaks with 
slowly varying centroid frequencies sometimes appear in the NS 
PDS for a few tens of minutes representing a timescale of 10* 
oscillations. We note that there is an intrinsic instrumental frag- 
mentation of observations that occurs on the same timescales and 
it is then assumed that the same QPO phenomenon often survive 
even longer (van der Klis, 2006). It would thus be neccesary to 
have an accretion flow with a neither low, nor high accretion rate 
that supports the considered torus-like configuration for a very 
long time. 

As quoted in Sect. 1, for the 3:2 epicyclic resonance the val- 
ues of the resonant eigenfrequencies for a non-geodesic flow 
are higher than those calculated for a nearly geodesic motion 
(Blaes et al., 2007; Straub & Sramkova, 2009). It has been 
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shown that in a certain case the difference can reach about 15%. 
This would change the non-rotating mass to a higher value, 
~ 1.2Mq. From Fig. 2 we can see that this value rather barely 
fits the modulation condition for the MIT EoS. It thus cannot 
be fully excluded that the model is compatible with observa- 
tions if the flow is fairly non-geodesic. However, needless to 
say that a serious treatment of this possibility will require in- 
vestigation of the related pressure effects on the disc structure 
in the Hartle-Thorne geometry, since the aforementioned studies 
only consider Schwarzschild (in a pseudo-Newtonian approxi- 
mation) and Kerr geometry. Moreover, they assume a constant 
specific angular momentum distribution within the disc, while 
there is evidence from numerical simulations of the evolution 
of accreting tori that real accretion flows tend to have rather 
near-Keplerian distributions (e.g. Hawley, 2000; De Villiers & 
Hawley, 2003). In this case, one expects the pressure corrections 
to be considerably smaller than those calculated for the marginal 
case of a constant angular momentum torus. However, there is 
no clear guarantee of this expectation and further investigation 
will be neccessary to resolve this issue. 

We can conclude that the resonance model for NS kHz QPOs 
should involve a combination of disc-oscillation modes that dif- 
fer from the geodesic radial and vertical epicyclic modes, or a 
modulation mechanism that differs from the Paczynski modula- 
tion. The results also suggest that the two modes together with 
the considered modulation may operate as long as a fairly non- 
geodesic accretion flow is assumed for strange- or some nuclear- 
matter EoS. 
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